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We propose an experimentally feasible scheme for topological interface engineering and show 
how it can be used for studies of dynamics of topologically nontrivial interfaces and perforation of 
defects and textures across such interfaces. The method makes use of the internal spin structure 
of the atoms together with locally applied control of interaction strengths to create many-particle 
states with highly complex topological properties. In particular, we consider a constructed coherent 
interface between topologically distinct phases of spinor Bose-Einstein condensates. 
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At the interface of two topologically distinct phases of 
a macroscopically coherent quantum system, the symme- 
try properties of the ground-state wavefunction change. 
Topological defects (e.g., vortices) cannot in general pen- 
etrate the interface unchanged. The boundary therefore 
has the property that defects must either terminate on 
the interface (typically as a point defect or monopole) or 
connect non-trivially to another object on the opposite 
side of the boundary. 

Interfaces between topologically distinct regions play 
an important role, e.g., in exotic superconductivity [1|, 



superfluid liquid 3 He A-B mixtures |2|-[4| and in Early- 
Universe cosmology. It has been proposed that a series 
of symmetry breakings lead to formation of domain walls 
and cosmic strings, which terminate on the boundaries 
between regions of different vacuum states [H|, 0] . Highly 
complex interface physics also arises from collisions be- 
tween branes [7] in string-theory brane-inflation jg] sce- 
narios. Superfluids have been discussed as candidates 
for experimentally accessible systems where analogues of 
cosmic topological defects may be studied [H, EJ For 
example, colliding liquid 3 He A-B interfaces have been 
proposed as analogues of string-theory branes [lOj . 

Here we show how atomic-physics laboratory tech- 
niques can be employed for engineering topologically 
non-trivial, coherent interface boundaries between spa- 
tially separated different ground-state manifolds which 
simultaneously exhibit different broken symmetries. We 
demonstrate non-trivial penetration of singular defects 
across a constructed stable interface between ferromag- 
netic (FM) and polar phases of a spin-1 BEC. We iden- 
tify the basic defect solutions crossing the interface and 
minimize their energies in order to characterize the de- 
fect core structures. We find examples of intriguing core 
deformations where a singular vortex terminates as an 
arch defect on the interface with the topological charge 
of a monopole, and where a coreless, non-singular vortex 
connects to a pair of singular half-quantum vortices. 



Our example demonstrates how the ultracold atom in- 
terface physics provides a novel medium for studies of sta- 
bility properties of field-theoretical solitons U The 
spin-1 BEC also already provides a possible system for 
dynamical investigation of phase transitions and defect 
production, e.g., of brane annihilation models. Moreover, 
the proposed method for interface engineering can exhibit 
especially rich phenomenology in spin-2 and spin-3 BECs 
and in strongly-correlated optical lattice systems. 

Atomic-physics technology provides tools for accu- 
rate detection methods for ultracold-atom systems on 
length and time scales that are difficult to achieve in 
laboratory systems of more traditional quantum fluids 
and solids. Advanced measurement techniques combined 
with the high degree of control over experimental pa- 
rameters make them suitable for quantum simulators of 
physical phenomena that are too complex even for nu- 
merical studies. This has attracted considerable interest, 
in particular in optical lattice systems, which can emu- 
late strongly-correlated condensed-matter models. The 
experimental development for using ultracold atoms as a 
laboratory testing ground for complex physical phenom- 
ena has been accelerated, e.g., by the observations of the 
Mo tt- insulator states of atoms lil-[l7|. the study of non- 
equilibrium defect formation in phase transitions in the 
Kibble-Zurek mechanism for both scalar [HI and spin-1 
BECs [19| , and in the preparation of artificial gauge-field 
potentials for multi- level atoms [2o| . 

Spinor BECs are condensates in which the spin degree 
of freedom is not frozen by magnetic trapping [2l| . They 
provide ideal models and emulators of complex broken 
symmetries due to their rich phenomenology of differ- 
ent phases 22l-[28| that support exotic defects and tex- 
tures [29l-l3l| . Spinor BECs have attracted recent experi- 
mental attention, e.g., in the studies of formation of spin 
textures (H, [Hj and in controlled preparation [HI and 
non-equilibrium formation [19| of vortices. 
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Here we will concentrate on a spin-1 BEC whose 
macroscopic wavefunction \£(r) may be written in terms 
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of the local density n(r) and a normalized spinor £(r) as, 



*(r) = V^)C(r) = VW) Co« , C f C = 1- (1) 

VC-(r)/ 

In the Gross-Pitaevskii mean-field description, the 
Hamiltonian density of the spin-1 BEC reads 



„_*| V ,|. + V(r)» + 2 B . + | rf i 
+ <?in^B-F^ + (; 2 n/(B-F) 2 



(2) 



(F) = ClFcxpCp i s the expectation value of the spin op- 
erator F, defined as a vector of spin-1 Pauli matrices. A 
weak external magnetic field leads to linear and quadratic 
Zeeman shifts described by the last two terms. The two 
interaction strengths are Co = 47rft 2 (2ao + (22) /3m and 
C2 = 4:7rh 2 (a2 — ao)/3m, respectively, where m is the 
atomic mass, and ao and CL2 are the s- wave scattering 
lengths corresponding to the two different values of the 
relative angular momentum of the colliding atom pair. 

The sign of C2 determines which phase is energetically 
favored by the interaction alone. For C2 < 0, as with 
87 Rb, minimization of the interaction energy favors the 
FM phase with the maximum spin magnitude, |(F)| = 1, 
in which case the broken symmetry of the ground-state 
manifold is defined by the rotations of the spin vector. 
The FM phase supports two topologically distinct classes 
of line defects 22, 231. The non-trivial vortices in each 
class are singly quantized singular line defects and non- 
singular coreless vortices, respectively (see Appendix for 



details), both of which have been observed 13, 132I. [3^. 

If instead C2 > 0, as with 23 Na, the interaction energy 
favors the polar phase, minimizing the spin magnitude, 
|(F) I = 0. The broken symmetry of the ground-state 
manifold is described by the unoriented nematic axis d 
(d = — d) and the condensate phase <\>. The polar phase 
therefore exhibits nematic order, analogously to liquid 
crystals and the superfluid liquid 3 He- A phase, and so 
supports both integer and half-quantum vortices. 

Scattering lengths in ultracold-atom systems are rou- 
tinely manipulated using magnetic Feshbach resonances. 
However, this is not possible in a spinor BEC, since the 
required strong magnetic field would freeze out the con- 
densate spin degree of freedom. It is possible to ma- 
nipulate scattering lengths also through optical (35| or 
microwave-induced Feshbach resonances [36] in which 
case the fields can be kept sufficiently weak in order not to 
destroy the spinor nature of the BEC. The Feshbach res- 
onance changes the scattering length by coupling the en- 
trance channel to a virtually populated bound state [36| . 
In particular, it is possible to tune the ratio ao/<22 be- 
tween the two scattering lengths. 

We suggest constructing an interface between topolog- 
ically distinct manifolds by spatially non-uniform engi- 
neering of the scattering lengths. In a spin-1 BEC this 



may be experimentally realized to prepare an interface 
between co-existing FM and polar phases. Using an op- 
tical Feshbach resonance, the spatial pattern correspond- 
ing to a sharp interface can be imposed using a holo- 
graphic mask. The spin-dependent interaction strength 
C2 is proportional to the difference between (12 and ao- 
Thus for small | C2 1 , as is the case with both 87 Rb and 
23 Na, only a small relative shift of ao versus CL2 is neces- 
sary to prepare the interface, and therefore the inelastic 
losses associated with optical Feshbach resonances [35[ 
can be kept small. 

A microwave field cannot similarly be focused. How- 
ever, using an optically induced level shift to tune the 
microwave transition off-resonant where no adjustment 
of the scattering length is desired, interactions may be 
manipulated in spatially well-defined regions to prepare 
a sharp interface boundary without the losses associated 
with the optical Feshbach resonance. 

In order to demonstrate the non-trivial nature of de- 
fect penetration across an interface between topologically 
distinct manifolds, we consider a harmonically trapped 
spin-1 BEC, where C2 abruptly changes sign at the cen- 
ter of the trap at z = 0. We have C2 > for z > 0, 
corresponding to the polar phase, and for z < 0, the 
BEC is in the FM phase with c 2 < 0. 

We first construct spinor solutions that approximate 
physical wavefunctions simultaneously in the two differ- 
ent manifolds and quickly relax to vortices connecting 
across the interface or terminating at the interface. Some 
illustrative examples of topologically allowed states are 
shown in Fig. [H The simplest connection can be iden- 
tified by considering a singly quantized vortex in both 
phases. Although such a vortex represents a very differ- 
ent topology in the two phases, it can be formed in both 
cases, e.g., by a 2tt winding of the condensate phase <p 
around the z axis. The two vortex solutions can be joined 
by changing the sign of either or £_ . By appropriate 
choice of parameters, doing so causes the spinor wave- 
function to adjust between the two manifolds by forcing 
|(F) I to switch from to 1 [Fig. QJa)], or else leads to a 
state which immediately relaxes to the desired configu- 
ration. Physically, such a sign-change in one of the two 
spinor components can be obtained by introducing a dark 
soliton plane (phase kink) in that component at z = 0, in 
which case the 7r-phase shift across the soliton is associ- 
ated with a vanishing density at the soliton core. As the 
density of the other components at the soliton core does 
not vanish, the BEC wavefunction continuously connects 
the two manifolds with the interface acquiring a width 
determined by the spin healing length £p = (871x2 n) -1 / 2 , 
the length scale over which | (F) | heals when locally per- 
turbed. 

A singular vortex with unit winding in the polar phase 
can also be written as a 2tt spin rotation about the z axis 
together with a 2tt rotation of the condensate phase </>. 
If we continue this solution to the FM side, by changing 
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FIG. 1: (color online). Schematic illustrations of possible vortex connections. The polar phase is above the interface and the 
FM below, (a) A singly quantized vortex in both phases, (b) A singly quantized vortex in the polar phase can connect to a 
coreless vortex or a Dirac monopole. These states are topologically equivalent. The monopole can be turned upside down to 
absorb the Dirac string in the polar phase, (c) A singly quantized spin- vortex terminates as a polar monopole. (d) A dipole can 
be constructed by joining the Dirac strings of a Dirac monopole and an anti-monopole [37j]. Placed on the interface, the dipole 
connects to two singly quantized vortices on the polar side, (e) The defect shown in (b) may be cut in half at the interface and 
the resulting vortices in the polar and FM regions may be moved apart if an additional dark-soliton plane is introduced in £o- 
More complicated vortex states may form as splitting of a singly quantized polar vortex into two half-quantum vortices (f), or 
by nucleation of half-quantum vortices that connect to coreless vortices that may exist together with monopoles (g). 



the sign of f_, we identify the resulting structure as an 
approximation of a coreless vortex [22], in which the spin 
profile quickly acquires a fountain-like texture. Hence we 
have constructed a solution where a polar, singly quan- 
tized vortex connects to a FM coreless vortex [Fig. [T]. 
We parametrize the spinor as 

/ -sin/3 \ 
c u+d v^cos/3 , (3) 

VZ \ ±e 2 ^sin/3 / 

where cp is the azimuthal coordinate, and f3 = 37r/4 gives 
an exact switch from polar to FM , with the negative sign 
in (_ used in the FM region. 

If the coreless vortex in Eq. (|3]) is continuously de- 
formed into a doubly quantized, singular vortex along 
the positive z axis terminating at the origin, the result- 
ing spin texture on the FM side forms a radial hedge- 
hog, (F) = r. This structure can be identified as the 
analogue of the Dirac monopole in quantum field the- 
ory, with the singular vortex line representing the Dirac 
string [3tJ [Fig. mb)]. The deformation is possible due 
to topological equivalence between the doubly quantized 
vortex and the vortex-free state. 

A singular spin- vortex in the FM phase can be made to 
terminate on a polar monopole [Fig.QJc)] as follows: The 
polar monopole is formed by two overlapping vortex lines 
of opposite circulation in £± perpendicular to a soliton 
plane in £n. The nematic axis then exhibits the radial 
hedgehog [29] d = f , which is the analogue of the t'Hooft- 
Polyakov monopole. Inserting a phase kink in £ + at the 
interface results in a structure on the FM side where (F) 
points radially away from the z axis, which we identify 
as the spin- vortex, 

/ Te -^sin<9\ 

where and cp denote the spherical angles, and the pos- 
itive sign is used in the FM part. 




FIG. 2: (color online). Core structure after minimizing the 
energy of Eq. J3}, corresponding to Fig. [TJb). (a) The mag- 
nitude of the spin (|(F)| = 1 is dark red). The polar vortex 
has split into two half-quantum vortices with FM cores with 
non- vanishing densities. White arrows show the spin vec- 
tor and indicate the coreless vortex in the FM part. Here 
c 2.0 x 10 4 ^j/ 3 , c 2 2.5 x 10 2 /kj/ 3 , Q 0.12cj and B = 0, 
where I = (h/muo) 1 ^ 2 is the oscillator length, (b) The nematic 
axis d (unoriented but shown as cones to emphasize winding) 
displays the characteristic tt winding as each half-quantum 
vortex is encircled. The two cores are joined by a disclination 
plane indicating the turn of d by tt. (02 = 1.0 x 10 4 ftwl 3 : 
n = 0.19a;.) 

So far we have analyzed the topological existence of 
defects perforating the FM-polar interface. In order to 
determine their energetic stability and the core structure 
we numerically minimize the energy of the constructed 
defect in a rotating frame F = E — Q(L Z ) by evolving 
the corresponding coupled spin-1 Gross-Pitaevskii equa- 
tions in imaginary time. Here Q denotes the frequency 
of rotation that is assumed to be around the z axis, 
(L z ) is the z-component of the angular momentum and 
E = j d 3 rV,(r). We assume a slightly cigar-shaped trap 
uj x = ujy = 2uo z = uj [46[ . 

Minimizing the the energy of C 1 ^ 01 results in the core 
deformation shown in Fig. [2j The state exhibits the core- 
less vortex on the FM side of the interface [Fig. (2a)] . The 
frequency of rotation determines the direction of the spin 
vector at the edge of the cloud. The core of the singly 







FIG. 3: (color online). Alice arch; when a singly quantized spin- vortex terminates on a polar monopole, the point defect 
deforms into an arch-shaped line defect, (a) An isosurface of the spin magnitude is shown in green. The spin magnitude rises 
to 1 (dark red) on the FM side of the interface (z < 0) and inside the Alice arch on the polar side. The singly quantized vortex 
with polar core remains in the FM phase. Grey cones indicate the spin vector, (b) The nematic axis (red cones) retains the 
hedgehog structure centered on the Alice arch (the deformed point defect), indicating that the topological charge is preserved, 
(c) The Alice arch is formed by deformation of the vortex cores in the f± spinor components. The vortex line in (red) splits 
at the position of the soliton plane. The upper part forms the Alice arch together with the vortex line in £_ above the interface. 
The vortex lines below the interface form the spin-vortex with polar core, (co = 2.0 x 10 4 /kj/ 3 , C2 = 5.0 x 10 2 /kj/ 3 , Q = and 
B = 0) 



quantized vortex in the polar part is deformed into two 
half-quantum vortices. 

The size of a singular defect core with vanishing den- 
sity is determined by £ n = (87rcon) -1 / 2 , the density heal- 
ing length. The singular polar vortex lowers its energy 
by spontaneously breaking axial symmetry, splitting into 
half-quantum vortices with FM cores of size £p. The de- 
formation is energetically favorable when the energy cost 
of the FM cores is smaller than the energy gained by re- 
moving the density depletion. We find the same splitting 
of the polar vortex for the defect with a singular vortex 
in both phases. 

A very intriguing core deformation results from mini- 
mizing the energy of a singular spin- vortex terminating 
on a polar monopole [Fig. HJc)]. The point defect re- 
quires the density to go to zero. For sufficiently large £p, 
the energy cost of the density depletion can be avoided 
by deforming the point defect into a semi-circular line 
defect with FM core whose ends attach to the interface. 
Fig. [3] shows the resulting arch-like defect. 

The arch is formed as a local deformation of the point 
defect, and retains its topological charge. Specifically, the 
radial hedgehog in the nematic axis d is preserved away 
from the arch. Single-valuedness of \£ then requires that 
d turn by tt on any closed loop through the arch, accom- 
panied by a 7r change in the phase <p. We thus identify the 
line defect as an arch-shaped half-quantum vortex, which 
we will call an Alice arch, as it is the interface analogue 
of the complete Alice ring [29]. Such ring-shaped defects 
are analogous to Alice rings in high energy physics [38| 
with a topological charge similar to magnetic "Cheshire" 
charge [39| We find that the deformation of a point de- 
fect to an arch is energetically favorable for C2 < 0.5cq. 



It is unstable towards drifting out of the cloud due to 
the density gradient, but could be stabilized by a weak 
pinning potential. 

Different defect structures penetrating the interface 
can be prepared experimentally by recognizing that the 
defects are composed of simple vortex lines, phase kink 
planes or, in more complex cases, by vortex rings in the 
three spinor components, each of which may be phase- 
imprinted using existing technologies [40-42]. Vortices 
may also nucleate due to rotation. We find that nu- 
cleation energetically favors defects consisting of a half- 
quantum vortex connecting to a coreless vortex, leading 
to states such as that illustrated in Fig. QJg). 

Here we have demonstrated topological interface en- 
gineering by studying examples of defect perforation 
across a FM-polar interface in a spin-1 BEC. The method 
can be extended to more complex broken symmetries in 
spin-2 [H, [U and spin-3 0, @ BECs that also 
support, e.g., non-Abelian defects [31(. Other particu- 
larly promising systems for topological interface phy sics 
are strongly-correlated atoms in optical lattices [15l4l7| 
exhibiting also quantum phase transitions and potential 
analogues of exotic superconductivity [1] in crystal lat- 
tices. 

Non-equilibrium defect production may be investi- 
gated in phase transitions in the presence of different 
broken symmetries [H, ^ . Defect formation from collid- 
ing interfaces can be employed as a model to simulate 
cosmological brane annihilation 0, H|. For instance, in a 
FM condensate, a slab of polar phase could be created, 
each interface being a 2D analogue of a D-brane. Re- 
moving the interaction shift causes the slab to collapse, 
bringing the interfaces closer until they meet and anni- 
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hilate, representing an annihilation of a brane-antibrane 
pair. In braneworld scenarios of cosmic inflation the an- 
nihilation may lead to defect production Q that could be 
directly observed in atomic BECs. A similar experiment 
has been performed with superfluid liquid 3 He A-B in- 
terfaces in which case, however, the detection of defects 
is difficult [loj]. 

We acknowledge discussions with D. J. Papoular and 
M. D. Lee and financial support from Lever hulme Trust. 



Appendix 

In this Appendix we present the basic defect configu- 
rations of spin-1 Bose-Einstein condensate (BEC) corre- 
sponding to the different broken symmetries of the ferro- 
magnetic (FM) and polar phases. Using these solutions 
we construct prototype spinor wavefunctions for the vor- 
tex connections across the interface between the polar- 
FM manifolds (such as those shown in Fig. 1 of the main 
text). The constructed prototype defect configurations 
are used as initial states for the numerical studies of the 
defect stability in the main text. 



Broken symmetries and vortices in the spin-1 
Bose-Einstein Condensate 



The wavefunction of a spin-1 BEC is written in terms 
of the density n(r) = |\I/(r)| 2 and a normalized three- 
component spinor in the basis of spin projection onto 
the z axis 



The spin-dependent interaction in spin-1 BEC gives 



C+(r) 



*(r) = V^WCW = Co(r) . (Al) 

VC-(r)/ 

The Hamiltonian density of a BEC of spin-1 atoms is [43| 



gm (B-F 



g 2 n 



((»•*)')■ 



(A2) 



where V(r) denotes the trapping potential, m is the 
atomic mass, F is the spin operator, given by a vector 
of spin-1 Pauli matrices, and (F) = CtFapCp- If a weak, 
external magnetic field is present, this will lead to linear 
and quadratic Zeeman shifts described by the last two 
terms. 

There are two interaction terms with strengths Co = 
47r/z 2 (2ao + <22)/3m and C2 = 47rft 2 (a2 — ao)/3ra, where a/ 
is the s-wave scattering length in the scattering channel 
with relative angular momentum /. The two correspond- 
ing healing lengths in the spin-1 BEC, £ n = (87rcon) -1 / 2 
and £p = (871x2 n) -1 / 2 , define the length scales over 
which n and |(F)| heal when locally perturbed. 



rise to two phases, the FM and polar [2 
l44j . Both phases exhibit distinct sets of possible defects. 
When a defect core is formed by a depletion of the den- 
sity, its size is determined by the density healing length 
£ n . A defect core in the polar (FM) phase may instead 
be filled with atoms in the FM (polar) phase. The size 
of the core is then determined by the spin healing length 
£,f [29(] . In the following we present the basic defect struc- 
tures in the FM and polar phases of a spin-1 BEC that 
are used in constructing the defect solutions crossing the 
coherent topological interface between the FM and polar 
phases. 

FM phase: If C2 < 0, the spin-dependent interaction 
will favor a state that maximizes the magnitude of the 
spin, i.e., |(F)| = 1. The general spinor may then be 
written as spin-rotations of £ = (1,0, 0) T , together with 
a condensate phase <\> 



V2e~ ia cos 2 f 
I = — — I sin^ 

/ ^2 I ^ e ia sin 2| 



(A3) 

where (a, /?, 7) are spin-space Euler angles, and 
7' = 7 — <j) absorbs the condensate phase. The broken 
symmetry of the ground-state manifold therefore cor- 
responds to the group of three-dimensional rotations 
SO (3). The spin vector is given by the Euler angles as 
(F) = (cos a sin /3, sin a sin /?, cos f3) . 

Topological stability of line defects can be character- 
ized by studying closed contours around the defect line 
and the mapping of these contours to the order param- 
eter space [45|. If the image of a closed loop encircling 
a line defect in the order parameter space can be con- 
tracted to a point, the defect is not topologically stable. 
SO (3) can be represented geometrically as S 3 (the unit 
sphere in four dimensions) with diametrically opposite 
points identified. The only closed contours that cannot 
be contracted to a point are then the ones connecting 
these diametrically opposite points, and we have only 
two distinct classes of vortices: singly quantized, singu- 
lar vortices that correspond to non-contractible loops and 
non-singular vortices representing contractible loops. All 
other vortices can be transformed to either one of these 
by local deformations of the order parameter. Examples 
of this first class are a singular line defect constructed as 
a winding of the condensate phase, 



cos 



2 p 



V2 



x/2 

sin f3 
V^sin 2 § 



(A4) 



where cp is the azimuthal angle, and the singular spin vor- 
tex constructed as a 2tt rotation of the spin vector 22|, Q 
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FIG. Al: Non-trivial vortices in the FM phase, (a) The sin- 
gular spin- vortex is formed as a radial disgyration of the spin 
vector (red cones) around the singular core, (b) Singular spin- 
vortex with cross-disgyration of the spin vector, (c) The core- 
less vortex is formed as a combined disgyration of the spin 
vector and a winding of the condensate phase, correspond- 
ing to spin rotations about the local spin vector (indicated 
by the orthogonal green and blue vectors). The core is non- 
singular and filled by the vortex- free spinor component. The 
direction of the spin vector at the edge depends on boundary 
conditions, (d) A different non-singular vortex configuration. 



[Fig. EUa) and (b)] 




(A5) 



A non-singular, coreless vortex is constructed as a com- 
bined rotation of the spin vector and the condensate 
phase [Fig. [Al{c) and (d)] 



V2 



cos 



2 P 



e i(fi sin/3 
V2e 2i * sin 2 | 



(A6) 



SO (3) does not strictly speaking support point defects. 
However, it is possible to form a spinor with a monopole 
structure of the spin vector (a radial hedgehog) as the 
termination of a doubly quantized vortex [37] [Fig. IA2j 
This is the analog of the Dirac monopole in quantum field 
theory. 



V2e~ 2 ^ cos 2 | 



sine 



V2t 



(A7) 



Polar phase: C2 > favors a state with |(F)| = 0. The 
general spinor is given by spin rotations of ( = (0, 1, 0) T 



FIG. A2: Dirac monopole; the analogue of the Dirac 
monopole is formed as the termination of a doubly quantized, 
singular vortex line (blue) resulting in a hedgehog configura- 
tion of the spin vector. 



together with a condensate phase 





C P = e^Z7(a,/3, 7 ) [ 1 





-e m sin/3 
V^cos/3 
e ia sin f3 



(A8) 

The unit vector d = (cos a sin /?, sin a sin /3, cos /3) defines 
the local direction of macroscopic condensate spin quan- 
tization. We may then write the spinor as 



icf) I d x + idy 
v/ V d x + id v 



(A9) 



Note, however, that C p (0> d) = £ p (</> + tt , — d). These two 
states must be identified in order to avoid double count- 
ing. The order parameter space, which might appear to 
be U(l) x5 2 , from the condensate phase and the rotations 
of d, actually becomes (U(l) x 6' 2 )/Z2, where Z2 denotes 
a discrete two-element group. The vector d is taken to be 
unoriented and defines the nematic axis, which together 
with (j) fully specifies the broken symmetry. 

A singly quantized vortex can again be constructed as 
a 2tt winding of the condensate phase 



c 1 




-e %a sin (3 
V2 cos 
e za sin /3 



(A10) 



The phase winding may be combined with a 2ir winding 
of the nematic axis by choosing a = ip 

— sin j3 



1 

71 



V2e ic ? cos/3 
sin (3 



(AH) 



2iip 



Because of the property that a ir change in <p can be ex- 
actly compensated by a tt rotation of d, it is possible to 
form vortices carrying half a quantum of circulation by 
letting both <fi and d wind by tt as the vortex is encir- 
cled [44[, leading to a disclination plane where d — ^ — d. 
If d is in the (x, 2/)-plane, a half-quantum vortex can be 
written 

-icp/2 \ / _l 
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(A12) 
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(b) 



1 ,L Vl> 



FIG. A3: Half-quantum vortex, (a) In the half-quantum vor- 



tex described by Eq. (|A12[) the nematic axis (red cones) winds 
by 7r as the core is encircled. This is accompanied by a ty 
winding of the condensate phase. The disclination plane is 
indicated by the dashed line, (b) Another possible tt winding 
of the nematic axis in a half-quantum vortex. 



In general, the axis about which d winds need not coin- 
cide with the vortex core. The ir winding of the nematic 
axis still allows us to identify the vortex [Fig. IA3I . 

The polar phase supports also point defects [29|. We 
can form a radial hedgehog of the nematic axis by choos- 
ing d = f = (sin#cos(/?, sin#sin(/?, cos#), which repre- 
sents a spherically symmetric monopole, 



1 

71 



suit 



cos 



(A13) 



suit 



Here the origin corresponds to a singular point defect. 
The two spinor wavefunction components £± form over- 
lapping, singly quantized vortex lines with opposite cir- 
culation. The vortex lines are oriented along the z axis, 
normal to a dark-soliton plane in the component (o- The 
structure of Eq. (|A13|) is the analogue of the t'Hooft- 
Polyakov monopole. 



Construction of prototype interface spinors 

We can construct an approximation of a desired de- 
fect configuration that crosses the FM-polar boundary by 
starting from a given FM (polar) defect state and insert- 
ing dark-soliton planes in some of the spinor wavefunc- 
tion components at the interface. For suitable choices of 
spinor parameters, the resulting configuration will then 
be sufficiently close to a local energetic minimum so that 
it quickly relaxes to the targeted defect structure. We 
will now explicitly construct the prototype spinors con- 
sidered in the main text. 

Singular vortex in both phases: To illustrate the 
method outlined above, we create a state with a singly 
quantized vortex on both sides of the interface [Fig. 1(a)] 
by starting from the FM singular vortex, Eq. (|A4|) , and 
inserting a soliton plane in (_ 



2 §_ 

2 



ei(p I V^cos- 
—= sin f3 



(A14) 



where the negative sign is used on the polar side of the in- 
terface. Note that for a general f3 the state created above 
the interface does not necessarily have |(F)| = 0, and 
therefore may not strictly be a polar state. However, the 
state created has the appropriate vortex structure and 
very quickly relaxes to the singly quantized polar vor- 
tex. The choice j3 = tt/2 yields a spinor which exactly 
changes from FM to polar. Analogous considerations ap- 
ply to other prototype spinors. 

Singular to coreless vortex: Inserting a soliton plane in 
(_ of the singular vortex C 1 in Eq. (|A11|) , results in the 
interface spinor 



(A15) 



where the negative sign is used on the FM side of the 
boundary. Choosing /3 = 7r/4 or /3 = 37r/4 yields |(F)| = 
1 on the FM side. The case f3 = 37r/4 corresponds to 
the coreless vortex ( cl [Eq. IA6j with fountain-like spin 
profile. 

Dirac monopole on the interface: A Dirac monopole 
may be placed on the interface by inserting a soliton 
plane into C D of Eq. (fA7|) 




V2e~ 2 ^ cos 2 f 



-up 



suit 



zby^si 



(A16) 



sin 



The spinor after the change of sign has a general structure 
similar to Eq. (|A11|) . The constructed spinor therefore 
approximates a singular vortex in the polar part, termi- 
nating on the Dirac monopole. Depending on which side 
of the monopole is placed on the polar side of the inter- 
face, the singular Dirac string either terminates on the 
monopole from the FM side, or becomes the singular po- 
lar vortex. The structure is equivalent to C X ^ C \ an d both 
spinors are represented by Fig. 1(b). 

Polar monopole on the interface: We consider a defect 
structure with overlapping, singly quantized vortex lines 
in £±, both oriented normal to the interface and of op- 
posite circulation, together with 7r-phase kinks in and 
Co- This can be parametrized as 



Asvfipm 




suit 



(A17) 



using the positive sign on the FM side. The resulting 
structure on the polar side represents a radial hedgehog 
of Eq. (jA13|) . The continuation of this on the FM side is 
similar to the singular spin- vortex £ sv in Eq. (|A5|) . Hence 
we have constructed an approximation to a spin- vortex 
on the FM side that terminates on the polar monopole 
on the polar side [Fig. 1(c)]. 

Terminating vortices: Vortices can also be made to ter- 
minate at the interface. In Eq. (|A15|h and in Fig. 1(b) of 
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the main text, a singular polar vortex perforates across 
the interface to a coreless FM vortex when we insert a 
7r-phase kink in The resulting defect can be cut in 
half while preserving the coherent interface with a con- 
tinuous order-parameter field across the interface. We 
achieve this by inserting an additional phase kink in (o 
that allows the vortices on either side of the interface to 
move apart 

/ -sin/3 \ 
C cut = -p \/2e^cos/3 , for z > (A18a) 
^ 2 \ e 2 ^sin/3 / 

_ / sin/3 \ 
C cut = -= V^e^cos/3 , for z < 0, (A18b) 
^ 2 \ e 2 ^sin/3 / 

where we may choose f3 = 37r/4 as in Eq. (|A15|) . This 
can yield the configuration in Fig. 1(e), where the singu- 
lar polar vortex and the FM coreless vortex are spatially 
separated and both terminate on the interface. Since the 
vortex lines in the individual spinor components termi- 
nate on the soliton planes, it is also possible to consider 
a state where a vortex exists only on one side of the in- 
terface, for instance, 

/ -sin/3 \ 
C pv = y/2e^cosP , for z > (A19a) 

^ 2 \ e 2 ^sin/3 / 

/ sin/3 \ 

C pv = — T V^cos/3 , for z < 0. (A19b) 
^ 2 V sin/3 / 
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